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Abstract 

The  ion  distribution  function  in  a  plasma  carrying  cxorrent 

(e.go  in  a  pinch  discharge)  is  not  Maxwellian,  and  this  fact  must 

be  taken  into  account  in  computing  the  thermonuclear  reaction  rateo 

In  particular,  it  is  the  (non-)Maxwellian  tail  of  the  distribution 

which  is  vital  to  this  computation,,   However,  the  conventional  method 

of  computing  the  distribution  is  legitimate  only  for  the  btilk  of  the 

distributiono   This  is  good  enough  for  the  computation  of  transport 

coefficients  but  not  for  the  computation  of  thermonuclear  outputo 

An  asymptotic  expression  is  derived  for  the  tail  of  the  ion  distribU" 

tion;  the  approximation  is  valid  when  the  ratio  of  electron  to  ion 

mass  Is  considered  to  be  small  and  the  electric  field  is  not  too 

large.   The  tail  end  of  the  distribution  is  foxmd  to  depend  exponent 

tially  on  the  value  of  the  electric  field,  and  depending  on  the 

numerical  values,  this  res^ilt  can  deviate  slightly  or  enormously  from 

the  Maxwell -averaged  reaction  rate.   However,  since  the  reaction  rate 

also  varies  very  rapidly  with  temperature,  an  estimation  of  plasma 

temperatiore  from  observed  values  of  neutron  emission  will  not  be 

greatly  in  error  even  if  the  computation  is  based  on  the  incorrect 

Maxwellian  distribution.      It   should  also  be  kept  in  mind  that  an 

experimental  temperature  deduced  in  this  way  will  be  the  maximum  that 

Is  attained  even  if  it  exists  only  in  a  very  small  region  of  the  plasmac 
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THERMONUCLEAR  REACTION  RATES  IN  A  DISCHARGE 
Harold  Grad 

1.   Introduction 

— ^^^— — — — —   • 

The  problem  of  evaluating  the  electron  distribution  (and  from  it 
the  conductivity)  of  a  plasma  which  carries  a  current  induced  by  an 
applied  electric  field  has  been  treated  by  Landshoff  [2]  using  the 
Chapman-Ens kog  approach  to  the  Boltzmann  equation  and  by  Cohen, 
Spitzer,  Routly,  and  Harm  [3]  using  a  Pokker-Planck  equation*   Exactly 
the  same  methods  could  be  used  to  obtain  the  ion  distribution  needed  to 
compute  thermonuclear  reaction  rates.   However,  it  is  easy  to  verify 
that  these  methods  yield  a  good  approximation  to  the  distribution 
fxinction  only  for  the  middle  range,  i.e.  for  ions  which  do  not  greatly 
exceed  the  mean  ion  speed.   A  symptom  of  this  defect  is  that  the  dis- 
tribution function  obtained  by  these  methods  is  not  necessarily  posi- 
tive.  One  can  expect  these  results  to  be  valid  up  to  a  maximum 
velocity  which  will  vary  as  1/E,   If  E  is  small  enough,  the  validity 
of  this  approximation  will  extend  out  to  the  range  where  thermonuclear 
reactions  occur;  however,  if,  as  is  frequently  the  case,  the  electrical 
discharge  is  used  to  heat  and  compress  the  plasma,  the  electric  field 
will  be  as  large  as  is  experimentally  feasible,  and  a  better  theory  is 
probably  necessary. 

The  theory  presented  here  is  also  a  linear  one,  applicable  to 


h  - 


small  values  of  the  electric  field,  but  it  is  a  different  lineariza- 
tion than  the  one  adopted  by  Landshoff  and  Spitzer  et  al.   Its  purpose 
is  to  obtain  a  uniform  approximation  to  the  distribution  function 
valid  for  large  velocity  as  well  as  small.   The  extension  of  these 
resxilts  to  larger  values  of  E  (i.e.  the  solution  of  the  nonlinear 
problem)  is  being  studied. 

A  simple  example  will  illustrate  some  of  these  possibilities  (but 
not  necessarily  in  the  way  in  which  they  occ\jr  in«the  problem  at  hand). 

Consider  the  ordinary  differential  equation  du/dx  =  Eu.   The  solution 

Ex 
is  u  =  u  e   .   The  linearization  u  =  u  +  Ev  leads  in  a  conventional 
o  o 

way  to  the  equation  dv/dx  =  Eu   with  solution  v  =  Eu  x.   This  is 
evidently  the  first  term  of  the  power  series  expansion  of  the  true 
solution  and  is  valid  only  for  small  values  of  Ex.   The  significant 
point  is  that,  even  in  a  linear  equation  in  which  a  parameter  E  occijrs 
linearly,  the  solution  is  not  necessarily  linear  in  E;  the  solution  ±a_ 
linear  if  E  multiplies  an  inhomogeneous  term.   In  the  conventional 
linearizations  of  the  Boltzmann  and  Pokker-Planck  equations  which  were 

just  mentioned,  not  only  is  the  quadratic  collision  term  linearized 

;9f  9f^°^ 

but  the  already  linear  term  E  •  — ^  is  altered  to  E  •  —  ^   where 

(o)  ^ 

f    is  the  zero-order  (Maxwellian)  approximation  to  the  distribution 

function  f.   We  propose  to  keep  the  term  E  •  ^y     \inaltered.   Also, 

in  the  collision  term  we  find  it  possible  to  introduce  a  linearization 

which  assxmes  only  that  there  is  a  small  departure  from  the  Maxwellian 

distribution  but  not  that  the  departure  is  a  small  fraction  of  f    « 

This  is  essential  in  approximating  to  the  tail  end  of  the  distribution 

where  f  can  be  strongly  non-Maxwellian  (but  still  small)  even  though  E 

is  small.   That  this  can  be  done  is  a  consequence  of  the  precise  form 

taken  by  the  collision  term. 
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In  this  context,  an  lllixminating  example  is  the  nonlinear  equa- 
tion du/dx  =  Eu  -  u^.   The  exact  solution  is  u  =  Eae^/d+ae^) . 
For  E  =  0  we  have  u  =  l/(x+b),  and  the  usual  linearization  is  clearly 
worthless  for  large  x.   However,  if  we  happen  to  know  or  guess  that  u 
approximates  E  closely  as  x  approaches  infinity,  the  linearization 
u  =  E  +  V  where  v  «  E  gives  the  correct  asymptotic  behavior  at 
infinity. 

In  this  report  we  have  chosen  the  Pokker-Planck  description 
rather  than  the  Boltzmann  description  of  a  plasma  for  mathematical 
simplicity.   Even  though  they  are  in  some  heuristic  sense  "almost  the 
same"  equation  (see  section  2),  for  our  purposes  the  Pokker-Planck 
equation  seems  to  have  distinct  mathematical  advantages.   Physically, 
the  motion  of  a  particle  (electron  or  ion)  in  a  plasma  is  governed  by 
the  superposition  of  many  small  simultaneous  random  forces  rather 
than,  as  in  ordinary  kinetic  theory,  by  a  time  sequence  of  distinct 
finite  tmpvilses.   A  particle  with  a  given  velocity  sees  a  directed 
frictlonal  force  tending  to  bring  it  to  rest  and  at  the  same  time  a 
random  dispersional  force  directed  "outward".   The  Maxwell  distribution 
results  from  a  balance  between  these  two  effects^   In  a  state  of  non- 
equllibrixjm  not  only  will  the  distribution  function  be  altered  but 
the  friction  and  dispersion  coefficients  will  also  be  perturbed. 

A  result  which  is,  at  least  at  first  glance,  surprising  is  the 
slight  extent  to  which  the  ion  distribution  is  altered  from  a  Maxwellian 
by  the  presence  of  an  electric  field.   Qualitatively,  this  may  be  seen 
as  follows.   Each  particle  sees  a  force  of  magnitude  +eE  and,  in 
addition,  a  frictional  force.   If  a  steady  state  is  postulated,  the 
mean  frictional  force  seen  by  the  whole  collection  of  electrons  (Ions) 


must  be  +eE  (-eE)o   This  mean  frlctional  force  Is,  of  course,  ascrib- 
able  only  to  the  alteration  of  the  friction  coefficient  concomitant 
with  the  non-Maxwell Ian  perturbation  of  the  distribution  function 
since  the  total  frlctional  force  vanishes  in  equilibrium.   In  this 
momentum  balance  we  need  consider  the  friction  on  an  electron  due  to 
the  ion  gas  only,  and  vice  versa.   For  the  electrons,  the  frlctional 
force  perturbation  is  a  definite  fionction  of  velocity,  and  the  differ- 
ence between  this  function  and  its  mean  value  (viz.  eE)  yields  the 
deviation  of  the  distribution  function  from  Maxwelliano   On  the  other 
hand,  because  of  the  great  disparity  between  the  two  masses,  the  ion 
velocity  distribution  is  almost  a  5-function  when  measured  on  the 
scale  of  the  mean  electron  speed.   Consequently,  the  friction  acting 
on  an  ion  due  to  the  pertxirbation  of  the  electron  distribution  function 
is,  to  a  certain  degree  of  approximation  the  same  for  all  ions.   Since 
the  mean  value  of  this  friction  must  be  eE,  and  the  friction  is  a 
constant  we  conclude  that  for  each  ion  the  applied  force,  -eE,  is 
exactly  cancelled  by  the  additional  friction  caused  by  the  perturba- 
tion of  the  electron  distribution.   Of  course,  this  Is  true  only  to 
a  certain  order  in  the  mass  ratio,, 

It  is  interesting  to  note  that  this  problem  can  be  formulated 

without  allowing  variation  in  time  only  in  the  linear  case,  for  small 

2 
E.   In  general  there  is  the  Ohralc  heating  term,  E<>  J  '-'  E  which  pre- 
cludes any  simple  stationary  solution.   There  is  a  special  case,  viz. 

with  a  time -dependent  electric  field  E  '^  t  '  -^ ,    which  yields  a  simi- 

2/3 
larity  solution  in  which  the  temperature  T  ^-^  t  '   and  the  distribution 

function  scales  as  '^'-^'^   f(4/  /T  ).   This  solution   has  not  yet  been 

developedo 
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2.   The  Fokker-Planck  Equation 

The  applicability  of  the  Boltzmann  equation  to  a  plasma  is  a 
considerably  more  subtle  affair  than  to  an  ordinary  gas.   The  reason 
is  the  presence  of  long  range  Coiilomb  forces.   Without  taking  into 
account  any  cancellation  of  forces,  one  concludes  that  the  force  due  to 
distant  charged  particles  overwhelms  that  due  to  neighboring  particles, 
This  is  in  fact  the  case  if  all  particles  have  the  same  sign  but  must 
be  modified  in  a  plasma  which  is  approximately  neutral.   The  first 
step  is  to  recognize  the  effect  of  the  distant  particles  as  being 
essentially  macroscopic  and  thus  to  be  included  in  Maxwell's  equations 
as  a  charge  density  and  current.   The  separation  between  "nearby"  and 
"distant"  is  not  arbitrary  and  must  be  taken  as  the  range  over  which 
a  significant  two-particle  correlation  exists.   This  correlation 
length  scale  is  the  Debye  length 


(2.1)     d 


2  ^    XkT 


2 


n  is  the  number  density  of  the  r'th  species  and  e   its  charge;  if 
the  ions  are  all  singly  charged, 

(2.2)     d2  =  ^  ,   n  =  Z  n^  . 

ne 

For  future  reference,  we  introduce  the  particle  diameter  (i.e.  mean 

distance  of  closest  approach) 

g2 
^2-3)      5   =   1271  XkT  ' 

We  use  mks  rationalized  units  except  where  otherwise  stated} 
X=   (3611  X  10^)"  is  the  permittivity  of  free  space. 
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the  mean  intermoleciilar  distance, 
(2.1^)  -K      =     4/3 

and  the  mean-free-path 

(2.5)  A  =  /^iog(d/a),  A^  =  ^  . 

Under  most   circumstances, 

(2.6)  5  «  X  «  d  «  y\_, 

and  a  particularly  useful  relation  to  remember  is  that  the  Debye 
length  is  essentially  the  geometric  mean  of  the  particle  diameter  and 
the  mean-free-path. 

Although  the  correlation  argument  satisfactorily  takes  care  of 
the  distant  (continuxim)  effects,  it  does  not  reduce  the  nearby  (colli- 
sional)  problem  to  conventional  kinetic  theory.   The  reason  is  that 
under  the  most  interesting  circumstances  the  Debye  length  is  large 
compared  to  the  intermolecxilar  separation  distance  and  this  fact  would 
seem  to  eliminate  conventional  binary  collision  analysis. 

To  resolve  this  difficulty  a  new  transition  length  somewhere 
between  5  and  X  is  introduced,  sayX.   For  distances  smaller  than  X, 
binary  collision  analysis  is  valid  and  the  Boltzmann  collision  term 
is  obtained.   In  the  remaining  range,  Z  <    x  <  d,  a  particle  sees  a 
superposition  of  many  small  forces.   They  are  small  in  the  sense  that 
each  one  alone  would  produce  a  small  deflection  in  the  particle's  path 
(x  »  5).   Following  Chandrasekhar  [ij  and  Spitzer  [3],  we  interpret 


Some  of  these  niimerical  factors  are  somewhat  arbitrary. 
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this  situation  stochastically  as  a  Markoff  process  in  which  each  of 
the  particles  in  the  range  ^  <  x  <  d  contributes  an  independent 
deflection  to  the  given  test  particle;  to  justify  the  Markoff  assump- 
tion there  is  also  assumed  to  be  no  correlation  over  times  longer 
than  the  mean  time  required  to  traverse  the  distance  d.   It  was  pointed 
out  by  Spitzer  that  this  stochastic  analysis  leads  to  the  same  result 
as  the  otherwise  unjustified  extension  of  binary  collision  analysis 
out  to  the  distance  d.   The  reason  is  that  the  Boltzmann  equation 
itself  represents  a  Markoff  process.   On  the  one  hand  (Boltzmann)  we 
compute  the  probability  distribution  of  binary  deflections  and  interpret 
them  as  independent  successive  impulses.   On  the  other  hand  (Fokker- 
Planck)  we  compute  the  same  probability  distribution  of  binary  deflec- 
tions and  interpret  them  as  simultaneous  independent  impulses.   Since 
the  mathematical  result  depends  only  on  the  probability  distribution 
and  not  on  the  interpretation,  the  Boltzmann  and  Pokker-Planck  results 
must  agree  (in  the  range  £<   x  <  d). 

Prom  this  discussion  we  see  that  the  Boltzmann  equation  is 
somewhat  more  accurate  than  the  Pokker-Planck  equation  in  that  it 
includes  large-angle  as  well  as  grazing  collisions;  the  quantity 
log(d/5)  compared  to  1  roughly  measures  the  importance  of  grazing 
collisions  as  compared  to  large  angle  collisions.   However,  in  choos- 
ing a  formulation  for  a  given  problem,  one  must  also  consider  the 
mathematical  feasibility  of  using  one  rather  than  the  other. 

The  Pokker-Planck  equation  is  derived  from  the  point  of  view  of 
the  basic  Markoff  process  in  [3]  by  Spitzer  et  al.  and  in  more  complete- 
ness in  [5]  by  Rosenbluth  et  al.   It  is  illuminating  to  indicate  a 
derivation  directly  from  the  Boltzmann  equation.   Por  simplicity. 
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consider   a   simple   gas. 


(^•7)       H  ^  ^r  H;  ^  ^  (V)  =  ^  J  ^  (^'^i  -  ^^i)  ff^  ^^^^^h 


(2.8)  K 


? ?5~T  • 


The  notation  is  the  conventional  onej  V  represents  the  relative 
velocity,  £,^   -  £,,    and  f  =  f(4'),  f[   =  f(4;[)  refer  to  the  velocities 
which  result  from  a  collision  between  particles  with  initial  velocities 
4,  ?^.   Specifically, 

f  ?«  =  C  +  a(a.V)  ,  a.V  =  V  cos  e 

(2.9)  < 

[  q  =  C^-  a(a.V) 

where  a  is  a  unit  vector  defined  by  the  cola  tit  ude  9  and  longitude  e 
taking  V  as  the  polar  axis. 

For  grazing  deflections,  the  differences  C  -  ^  and  4^  -  ^  are 
small  and  we  may  expand  the  difference  f'f'  -  ff_  in  a  Taylor  series 
using 

(f'    =  f  +  V  cos  9  (a.V)f  +  ^  V^cos^9(a.V)^f  +  ».o 

(2.10)  <  , 

\f[   =  f ^  -  V  cos  9(a.'^)f^  +  ^  V'^cos'^9(a.\ir)'^f^  +  ... 


Carrying  out   the   e    integration,    we   get 

r(f'f^   -  ff^)de   =  Tt  cos^9  ((1  +  I  sin^G)(V-\7     -  V«  V,)ff, 

(2.11)  2  ?  1 

+    (1   -  I   sin''e)(VV-  V.7L)'^ff^   +...{-     . 

Here  V=  ^/^C  and  V^  =  9/d^^.      The  9  integration  diverges.   This 
difficulty  can  be  avoided  either  by  using  the  correct  Debye  shielding 
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potential  instead  of  the  Coiilomb  potential  or  by  the  conventional 
expedient  of  introducing  a  cutoff  at  an  angle  9  which  corresponds  to 
the  deflection  produced  by  a  particle  approaching  with  mean  relative 
velocity  at  a  separation  d<,   The  result  is 
(2.12)     cos  9^  =  5/d. 

Two  integrals  arise, 
9 


/ 


cos^9  .5l2_®  d9  =  -log  I  cos  9  | 
cos-^9  ° 


J 


cos^9  sin^9  ^^^  d9  =  -log]  cos  9^  1  =  ^   „ 


cos-'^9 


The  1/2  is  dropped  compared  to  the  dominant  term  log(d/5)  since 
equivalent  terms  have  already  been  dropped  by  the  elimination  of 
large  angle  deflections.   Moreover,  the  third  and  higher  derivatives 
in  the  Taylor  series  employed  above  can  be  ignored  only  by  using  this 
rule.   After  many  integrations  by  parts,  we  obtain  the  result 

(2.13)   H  +  ^T.  ^  -^  5I-  (V)  =  i-  ^Ag  ^KJ^  "  rf-  (^r.^ ) 

^t    r  ^x^   t?Cj,   r     2  o'C^<?^g   rs     ^E,_^        r 


where 


(2.11;) 


J    /  V  O    "V  V 


We  recognize  a^  to  be  the  friction  coefficient  and  b   the  dispersion 
coefficient.   A  particle  which  is  know  to  have  the  velocity  C  at  time 

r 


R.  Liboff,  report  to  appear. 
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zero  "diffuses"  in  the  Interval  dt  into  a  normal  distribution  with 

mean  value  £,     +  a  dt  and  second  moments  b  dt. 
r    r  rs 

For  a  mixture  of  electrons  and  several  sorts  of  ions,  we  quote 
the  results 

Where 

(2.16)  ai(4)  =  L^  Z  ^^^-V-  Ai(€) 

(2.17)  b^^  (4)  =  L^Tb^^(C) 

(2.18)  L^  =  -4?  l°6(f ) 

(2.19)  A^(C)  =  r^  f^(';)d>(         (Vp  =  ^j.  -  ^^) 

,        /^  V^5    -  V  V 

(2.20)  B^g(?)  =  J  ^°  3  ^  ^  f J(7)dy 

These  results   are    identical   to   those   given  by  Rosenbluth  et   al. 
in   [5]. 
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3o   Properties  of  the  Friction  and  Dispersion  Coefficients 

For  reference  we  list  a  number  of  useful  identities  satisfied 
by  the  various  coefficients: 


(3.1) 


J  aJ(^)  fJ(C)  d^  =  0 


(3.2)     21  ^^  J  aj(^)  f^(?)d4  =  0 


(3.3) 


I  Ik^   A^(4)  +  bJ^(4)]  f^(€)d4  =  0 


(3.1|)     Z:  ^i  J  t24paj(?)  +  bj^  (C)]  f^(4)d4  =  0 


i 

5B^        . 
(3.5)     _Zi=2AJ 

^^s 
The  first  two  identities  are  useful  in  deriving  the  momentum  conserva- 
tion equations  and  the  next  two  are  similarly  related  to  energy  con- 
servation.  The  final  identity  will  be  foxmd  usefxil  in  manipxilations 
to  follow.   The  proofs  of  these  Identities  are  quite  direct. 

We  shall  find  it  useful  to  expand  the  indicated  derivatives  in 
the  collision  term,  J.,  in  (2.1 5), 


(3.6) 

jl 

=  i^^3l^-Jii-ft^^ 

where 
(3.7) 

4 

-'k-r-^^r'y4 

(3.8) 

< 

-^|;^=.vz:^^^ 
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The  coefficients  P  and  Q  are  now  expressed  In  terms  of  A  alone » 
A  significant  point  is  that  the  P   stmmiation  contains  only  terms 
with  j  7^  i. 

The  various  coefficients  can  be  explicitly  evaluated  for  a 
Maxwelllan  f  (cf.  [l]),  viz. 


(3.9) 


.(0) 


n ^-4.V2RT 


Tf^^ 


R  =  k/m 


(aiRT)- 

(the  superscript  i  has  been  omitted).  By  Introducing  V  as  the 
variable  of  Integration  instead  of  >?  and  introducing  polar  coordinates 
with  the  vector  ?  as  the  polar  axis  we  obtain 


P  (-^) 


(3.10) 


where 


(3.11)         < 


rif^°^y))d7    = 

rvf^°^r^)drj     =  n  )/Rf     G(-^) 


G( 


„=/i    [,x.  1)   |\-<l/2)7%  .   e-<l/2)'='| 


By  differentiation  we  then  find 
(3.12)     a(o)=^  ^P'(X) 


(3.13)    B. 


(o) 


/If  L 


rs    r  s 


G  (x)  +  ^^^0  (x) 

X 
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where  we  have  used  the  abbreviation 
(3.11;)       ^r  =  7=  • 
For  small  £,   we  have, 

(3.16,   B<:)=/|-^[a„(|-^^.^.'^....) 


>^RT 


/   2  _^  1   2  ^ 
'^s^-  13-"  3F^  -^ 


•'] 


and  for  large   5j 


r   2  2 

,v                  lx5-xx              x5        -3xx 
(3.18)  B^°)=-;g^^  rs^^rs ^^  ^        ^  ^ 


/RT 


X'  X' 


the  terms  neglected  are  in  this  case  exponentially  small. 

In  the  presence  of  an  electric  field  f  is  no  longer  Maxwellian. 
For  a  perturbation  which  is  linear  in  E  we  have 

(3.19)    f(?)  =  f^°Uv    +  E.^  0i^   ?2). 

By  the  assTomption  of  linearity,  ^  is  a  scalar  function  of  the  vector  £, 
alone  (not  of  E)  and  is  therefore  a  function  of  only  the  magnitude  of 
C«   For  the  perturbations,  defined  by 

'^•^°'  ^      b'°)  .  b'1) 

rs       rs 
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we  have 

(3.21)  A^^^  =  Eg  (II^  0{^ri^)AYl 

and  2 

(3.22)  B^^^  =  E^  J  ""  ^^^3  ''^''"  r|^  J2^(|  n^)d>^ 

We  are  particularly  Interested  in  the  values  of  A    and  B    for 
small  £,,      To  this  end  we  expand 

(3.23)  =  0{^  v^)  +  v.C  J^'(|  v^)  +  i  4^0' (i  v2)+  |(v.c)V'(|v2) 

+  l^^v.^  /   (|v^)  +  ^(v-c)^  /'(|v2)  +  ... 

obtaining,  after  some  manipulation 

(3.21^)     A^l)  =^n  0^E^   -  |u  0^(E^52  ^  2E.?  4^)  +  ... 

(3.2^)     B^l^   =  ^0^{E.^  6^^   -H  E^^g  +  Eg^) 


where  we  have  set 


0^      =      0(0) 
(3.26)   <        .«> 

1^0^   =  J   0(x)dx 

o 
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1;,   Reduction  to  Spherical  Coordinates 

In  our  problem  the  external  force  terra,  X   In  equation  (2,15), 
Is  the  electric  field. 


(l|.l) 


T  ^r' 


We  can  use  spherical  or  cylindrical  coordinates;  we  choose  spherical 
coordinates  (r,0)  with  the  direction  of  E  as  the  polar  axis. 


(1^.2) 


C  =  Er  cos  0. 


In  order  to  transform  the  collision  term  in  the  form  (3.6)  we  note 
that 


M  -5^      A  -   ^^   _,p^.2^.1   <?".  cot  9  _^ 


also 


=r  <?£ 


d^ 


C^ 


(i^.14-)  < 


r-s  9K,^e^^ 


,2  Jl 
9v' 


-r.-fc  =  ^(-«^-^^  f.^ 


L  (^r^s-  Vr-  2E.5  ^  )  ^^  =  -2E  sin  9(^  -  i  ^ 


Next  we  set 

r 


?^^, 


b   =  b. 


rs    1   ^2 


+  b^(5 


^«€« 


C^5. 


2'  rs 


-^  )  -b3(E^€3.E^€^-2E.^^  ) 


(1+.5)  <   Pp  =  Pl^,.  +  P2I 


Qj.  =  qi(r)?r  +  q2(r)E.^  4^  +  q3(r)E^ 
=■  19  - 


and  compute 

(1,.6)     ^  =  3q^  .  r  ^  +  E  cos  0(Upq2  -  r^  ^     +  ^  ), 

Note  that  b    is  separated  Into  three  orthogonal  components  in  the 
sense  that  the  product  of  £  S;,  /^   and  5   -  £,£,/£,      is  zero,  etc. 
Dropping  the  time  and  space  derivatives  in  the  Pokker-Planck 
equation,  we  obtain 

+  [il,2  +  rp^  +  Eoos9(p2-l)]|£-^ 
(1|.7) 

+[  Ibg  cot  9  ^-  E  sm  e(b3  -  pg  4.  |)]  i  II    - 

-[  311  *  "^  +  E  CO.  9(1^,2  -^  r2  ^  .  ^)]f  =  0, 
the  superscript  i  has  been  omitted. 
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5.   Solution  for  Thermal  Speeds 

In  order  to  carry  out  the  asymptotic  evaluation  of  f  for  large  4 
in  the  next  section,  it  is  necessary  to  first  discuss  the  behavior 
of  the  solution  for  moderate  values  of  4.   In  this  range  we  can  use 
the  conventional  linearization  (3.19).   It  will  not  be  necessary  to 
solve  this  problem,  only  to  estimate  the  order  of  magnitude  of  the 
solution.   Substituting 

(5.1)     f  =  f ^°^  +  f^^^  =  f^°^  +  E  cose  r  0 

into  equation  (I;.?)  and  linearizing  with  respect  to  E  we  obtain 

l^(o)  Zfi^^  1^(0)  A-^ffiil^  (l^(o)  ^^^(o))£f^ 


.(1)   .   r.^     <9qi°^ 


(5.2) 


+  ib(°)cot  9  i  ^  -(3q^°^  +  r  1^  ) 

+  ^b^   cot  9  ^  ^e     ^^^1     ^  o'r   ^ 

1  ^(1)/  r^     1  X    /lv(l)  +  „„(l)^^  +  r-  ^„c  Qf©   T>  "1^ 


2  -1   R.^. 


;4l)..^^.Kcos9(Ur4o)_2^,^,- 


f(o) 


The  inhomogeneous  terms  have  been  collected  on  the  right  hand  side. 

It  is  easy  to  verify  that  E  cos  ©  cancels  out  as  a  common  factor 

-1   p 
leaving  an  ordinary  differential  equation  for  0{-^   r  ).   We  wish  to 

estimate  the  magnitude  of  the  various  Inhomogeneous  terms  considering 

that  m'/m"^  =  m/M  is  small.   The  given  equation  (5.2)  Is  supposed  to 

be  that  for  the  ions.   Consequently,  we  are  Interested  in  a  range  of 
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velocities  r  comparable  to  ^kT/M  .   The  dominant  term  on  the  right 
(this  arises  from  the  original  term  E*  af/^C)  is 

E  cos  9  I  ^   . 

Now,  following  Spitzer,  we  can  evaluate  the  parameter  p-  which  sub- 
tracts from  this  dominant  term  by  referring  to  the  conservation  of 
momentum.   It  is  easiest  to  use  the  form  (2,1^).   Multiplying  through 
the  ion  equation  by  4p  and  integrating  over  velocity  yields 


i 


(g  E^  +  a^(?))  f''(5)dC  =  0. 


Using  (2.16)  and  the  identity  (3.1)  for  a"^  we  are  left  with 


Ja;(4)  f-^i 


I  E^  n-"  +  L-*-  ^  I  A-ja    f(C)d4  =  0. 


Linearizing  A~  and  f   in  the  integrand,  we  find  that  the  A   f 

term  is  small  of  order  m/M  compared  to  the  A   'f^®^  term  and,  using 

the  first  term  in  the  expansion  (3«2l4.)  (the  whole  ion  distribution 

is  at  the  very  low  velocity  end  of  the  electron  distribution)  we  obtain 


Now,  evaluating  P  from  (3o7)  and  (3»2l4.)  and  from  this  the  coefficient 


P2  (cf.  ik'^)),   we  find 


r        m    r 


(5.U)      p^.-^^-L^S  =  I 
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We  see  then  that  what  we  would  expect  to  be  the  dominant  terra  in  the 
equation  drops  out,  the  reason  being  the  one  described  in  the  intro- 
duction, namely,  the  cancellation  of  the  applied  E  by  the  pertxirba- 
tion  of  the  ion-electron  friction  coefficient.   It  is  very  easy  to 
verify  that  the  other  inhomogeneous  terms  in  equation  (5o2)  are  also 
small,  on  the  order  of  m/M;  the  coefficients  b,   ,  etc,  are  obtained 
by  using  the  low  velocity  approximation  to  the  electron  coefficients 
derived  in  section  3.   It  would  be  possible  to  solve  (5«2)  just  as 
in  [3]  but  this  will  not  be  necessary.   As  it  turns  out,  only  the 
coefficient  0Z   of  the  pertxorbed  electron  distribution  is  needed  (and 
this  has  just  been  evaluated  In  (5.3)).   The  perturbation  of  the  ion 
distribution  in  the  thermal  range  can  be  ignored  to  low  order  in  an 
expansion  in  m/M„ 
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6.   Asymptotic  Expansion  of  the  Tall  of  the  Distribution 

We  Introduce  as  an  expansion  parameter  the  mass  ratio. 


e^ 


i '    M 


R  T 

For  deuterons  M/m  '-'  3600,  CVc"*"  ~  60,  and  e  ^  l/l;.   Actually  we  have 
two  small  parameters,  E  and  e,  but  they  will  eventually  be  related. 
We  are  interested  in  a  range  of  ion  velocities  large  compared  to  C 
and  small  compared  to  C~,   In  this  intermediate  range  the  Maxwellian 
friction  coefficient  a"^  has  the  approximate  form  (from  (2.16),  (3.1^), 

(3.17)) 

a  =  a   +  a 


=  2L^A+   +  \  l\- 
m 


+ 

R"T 


^<%^i/l^'-' 


(6.2) 

where 
(6.3) 


and  the  absolute  magnitude  of  the  vector  a  has  been  taken.   The 
minimum  value  of  a   is  found  at 

(6ol4-)     ex  =  (7201)^/^--  2.5  . 

At  this  point  the  value  of  the  friction  coefficient  is  reduced  by  a 

2  H 

factor  on  the  order  of  e   times  the  peak  values  attained  near  ^  '^  C 

and  £,  r^C~, 
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It  is  not  clear,  a  priori,  whether  we  are  Interested  in  the 
range  x  --v  l/e  or  x  /v  l/e   (friction  Ay  e   or  e ) ,   Presumably  the 
correct  scaling  depends  on  which  range  of  velocities  produces  the 
thermonuclear  reactions  in  a  given  experiment;  however,  since  there  is 
only  a  factor  of  four  at  stake,  mathematical  expedience  would  seem  to 
be  the  predominant  criterion. 

Several  Important  introductory  remarks  will  be  made  now.   First, 
for  small  E,  even  though  the  tail  end  of  the  distribution  may  be 
extremely  non-Maxwellian  it  will  still  die  down  exponentially.   Now, 
the  friction  and  dispersion  coefficients  are  given  as  integrals  over 
the  whole  distribution.   They  will  therefore  be  determined  (even  the 
values  at  high  velocities)  by  the  middle  range  of  the  distribution 
for  which  the  conventional  linearization  (3.19)  is  valid.   The  use  of 
this  linearization  in  computing  the  perturbations  (3o2I^.)  and  (3o25) 
is  in  this  way  justified.   For  those  coefficients  which  are  integrated 
over  the  electron  distribution,  we  can  use  Spitzer's  determination  of 
this  distribution.   For  the  ion  coefficients  we  can  take  the  iinperturbed 
Maxwellian  values  to  the  order  which  will  be  kept  in  the  expansion 
(cf .  section  5) « 

The  reason  why  the  present  linearization  gives  correct  results 
even  if  the  perturbation  is  large  compared  to  f   ,  provided  only 
that  it  is  small  in  an  absolute  sense,  is  that  instead  of  expanding 
a  quadratic  functional 

Q(f,f)  =  Q(f^°\  f^°^)  +  Q(f^°^  f^^^)  +  Q(f^^\f^°^)+Q(f^^\f^^^ 

and  then  dropping  Q(f  ^ -"-^f  ^-"-^ )  compared  to  Q(f  ^°^  ,f  ^■'"^ ) ,  we  find  it 
possible  to  expand 
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Q(f,f)  =  Q(f,f^°^)  +  Q(f,f^^^) 

where  the  f    has  been  found  independently  and  is  not  an  unknown. 

Adopting  these  approximations,  we  obtain  the  following  equation 
valid  in  the  range  C  «  £,  «  c"   after  factoring  out  the  Maxwellian, 

(6.5)     f  =  f^°^g. 

It  is  easy  to  see  by  inspection  that  taking  E  ^  e   and  r  -^  l/e 
yields  a  solution  for  g  which  differs  by  only  a  small  amoiont  from 
unity;  the  conventional  linearization  would  therefore  apply  out  to 
this  range;  (we  must  keep  in  mind  that  0-.    has  order  e  ).   On  the 
other  hand,  taking  r  ^l/e    ,  we  see  that  the  terms  in  9g/<?r  and  g 
are  dominant  and  we  obtain  the  asymptotic  solution 

(6.7)     g  =  exp  J  -  Hji  /IpF2^^|^ f^^^         E  cos  9  i  . 

1    ^      ne-'log(d/5)  V2R  T/  I 

The  exponent  has  order  \inity  when  the  energy  gained  by  a 

particle  accelerating  under  E  for  a  mean-free-path  is  comparable  to 

kT. 

For  E  '^  e  and  r  ^  l/e      the  exponent  in  (6.7)  has  order  l/e^. 
This  should  be  compared  with  the  exponent  in  the  Maxwellian  which 
wooild  then  be  of  order  l/e^.   Clearly,  for  e  small,  the  perturbation 
is  a  very  large  factor  applied  to  the  Maxwellian,  but  the  Maxwellian 
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is  extremely  small  to  an  even  higher  order.   The  asymptotic  expansion 
for  g  has  been  carried  to  only  one  significant  term;  the  factor  g  is 
not  acciarate  to  this  order  but  its  logarithm  is.   Fortunately  this  is 
sufficient  to  obtain  interesting  information. 

From  the  form  taken  by  g  In  (6,7)  it  is  quite  clear  how  to  pro- 
ceed to  get  a  uniform  approximation  to  f  for  thermal  as  well  as  large 
velocities.   The  linear  perturbation  f  =  f  ^^  +  Er  cos  ©  0{jv    )  has 
been  seen  to  be  valid  up  to  r  -^l/e.   But,  so  long  as  the  perturbation 
is  small,  this  is  entirely  equivalent  to  f  =  f   exp  \e   cos  ©\}>(r)| 
where  r^  =   f   \Jf,   This  exponential  form  is  therefore  valid  for  the 
whole  range  of  r  (up  to  r  /^  l/e  ),   If  it  is  substituted  into  the 
Pokker-Planck  equation  and  the  result  linearized  with  respect  to  E, 
an  ordinary  differential  equation  for  )|f(r)  will  result,  similar  to 
the  one  solved  by  Spltzer  [3],   This  analysis  will  not  be  carried  out 
here  since  the  asymptotic  formula  is  probably  svif f Iciently  accurate 
for  the  present  purpose.   It  should  be  noted  that  although  this  asymp- 
totic expansion  can  be  continued,  computing  additional  higher  order 
terms,  there  will  remain  an  undetermined  constant  of  integration 
(multiplicative  factor  on  g).   This  constant  can  only  be  determined 
by  matching  the  asymptotic  solution  to  the  low  velocity  result,  i.e. 
by  carrying  out  the  uniformly  valid  procedure  just  mentioned. 

It  should  be  relatively  easy  to  detect  the  presence  of  the 
phenomenon  described  by  equation  (6,7)  because  of  the  remarkable 
fact  that  the  deviation  from  Maxwellian  is  opposite  to  what  would  be 
expected  from  the  direction  of  E.   This  is  not  surprising  when  we 
recall  that  the  presumably  dominant  term  in  the  Fokker-Planck  equation, 

Of 

E»  -T^   ,  was  cancelled  by  the  perturbation  of  the  electron  friction 
coefficient. 
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7.   Thermonuclear  Output 

We  can  vn?ite  the  asymptotic  expression  for  g,  equation  (6.7), 
as 

(7.1)      g  =  exp 
where 


[-a(E.4)5^j 


(7.2) 


a  =  ^  ^3/2  _^ 


eV/2 


ne-^log(d/5)  /2kf" 


To  dominant  (i.e.  exponential)  order,  the  D-D  cross-section  takes 
the  form 

(7.3)      o-  ~exp  |ti(lMV^)  j  . 

To  dominant  order,  the  average  cross-section,  i.e.  the  thermonuclear 
output,  is  given  by 


K 


lMr2 


(7.1;) 


=  J  exp  k(i  MV2)  -  2__  _  e,(E.4)^^-  %^  -  d(E.4-L)cJ 
=  r  exp  L(|mv2)-  jj^(V^  +  W^)-  ^EW(3v2  +  W^ )  1  d^d^^^ 


d€d?, 


where 

(7.5)      W  =  C^  +  5, 


V  =  ?,  -  ^. 


To  the  same  dominant  order,  the  value  of  v  can  be  approximated  by 
maximizing  the  exponent  with  respect  to  V  and  W.   This  occurs  at 
values  V^ ,  W-  given  by 

r 


(7.6) 


W^.-lf  aEv2 


'   u'(iMv2)  =  -1-  -  9  kT  ^  V^ 
H.  i^v^;   2kT  %  ^^  ^^  ^o 
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(V,  is  given  implicitly  by  the  second  relation).   Here  V   is  the 
maximizing  value  for  the  Maxwellian  distribution,  i.e. 

(7.7)  ^(^^)  =  ^  . 

Inserting   (7,6)    into    (7 •14-)    and   treating  E   as   a  small   pertvtrba- 
tion,    we   obtain 

(7.8)  v  =  v„e.p|^f  aVv^} 

where  v   is  the  Maxwellian  reaction  rate, 
o 

For  temperatures  which  are  not  too  high  we  have  (e.g.  see 
Thompson,  [l\.]  ) 

kT     T  /^ 
on  the  right,  the  temper  attire  T  (=kT)  is  in  kev.   Rewriting  (7o8)  in 
these  units,  we  obtain 

where  n  is  in  ions/cm  and  E  is  in  volts/cm.   For  example,  taking 
n  =  10-^\  T  =  0.5  (5  X  10^  °C.)  for  which  log(d/5)  ^  ll^,    we  find  that 
the  correction  becomes  appreciable  for  E  ~'  1  volt/cm.   Any  larger 
value  of  E  will  produce  an  enormous  increase  in  thermonuclear  output. 
The  most  striking  feature  of  this  formula  is  the  extremely  sensitive 
dependence  on  E,  T,  and  n.   We  note,  however,  that  the  perturbation 
from  V   to  v..  is  small  even  when  v  as  computed  from  (7.8)  is  much 
larger  than  v  .   This  means  that  use  of  the  Maxwellian  formula  to 
estimate  T  from  an  experimentally  observed  value  of  v  will  not  produce 
a  large  error.   One  must  be  caj?eful,  however,  to  interpret  such  a 
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temperature  as  the  maximum  which  exists  locally  anywhere  In  the 
plasma. 
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